Entanglement and nonclassicality of photon-added two-mode squeezed thermal state 
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We introduce a kind of entangled state — photon-addition two-mode squeezed thermal state (TM- 
STS) by adding photons to each mode of the TMSTS. Using the P-representation of thermal state, 
the compact expression of the normalization factor is derived, a Jacobi polynomial. The nonclassi- 
cality is investigated by exploring especially the negativity of Wigner function. The entanglement 
is discussed by using Shchukin-Vogel criteria. It is shown that the photon-addtion to the TMSTS 
may be more effective for the entanglement enhancement than the photon-subtraction from the 
TMSTS. In addition, the quantum teleportation is also examined, which shows that symmetrical 
photon-added TMSTS may be more useful for quantum teleportation than the non-symmetric case. 

PACS number(s): 42.50.Dv, 03.65.Wj, 03.67.Mn 



I. INTRODUCTION 

Quantum entanglement with continuous- variable is an 
essential resource in quantum information processing 
[l|, such as teleportation, dense coding, and quantum 
cloning. In a quantum optics laboratory, a Gaussian two- 
mode squeezed vacuum state is ofen used as entangled 
resource, which cannot be distilled only by Gaussian lo- 
cal operators and classical communications due to the 
limitation from the no-go theorem [2|-Q- To satisfy the 
requirement of quantum information protocols for long- 
distance communication, there have been suggestions and 
realizations for engineering the quantum state, which are 
plausible ways to conditionally manipulate a nonclassical 
state of an optical field by subtracting or adding photons 
from/to a Gaussian field |5l-[ll]. Actually, the photon 
addition and subtraction have been successfully demon- 
strated experimentally for probing quantum commuta- 
tion rules by Parigi et al. |12| . 

In order to increase quantum entanglement, two- mode 
photon-subtraction squeezed vacuum states (TPSSV) 
have received more attention from both experimental- 
ists and theoreticians |5l, fiol, ll3l - [22] . Olivares et al. [i3| 
considered the photon subtraction using on-off photo de- 
tectors and showed improvement of quantum telepor- 
tation, depending on the various parameters involved. 
Kitagawa et al. [14], on the other hand, investigated 
the degree of entanglement for the TPSSV by using an 
on-off photon detector. Using an operation with single 
photon counts, Ourjoumtsev et al. [ill, [l6| demonstrated 
experimentally that entanglement between Gaussian en- 
tangled states can be increased by subtracting only one 
photon from two- mode squeezed vacuum states. In addi- 
tion, Lee et. al [22] proposed a coherent superposition of 
photon subtraction and addition to enhance quantum en- 
tanglement of two-mode Gaussian sate. It is shown that, 
especially for the small-squeezing regime, the effects of 
coherent operation are more prominent than those of the 
mere photon subtraction and the photon addition. 

Recently, we proposed the any photon-added squeezed 



thermal state theoretically, and investigated its non- 
classicality by exploring the sub-Poissonian and nega- 
tive Wigner function (WF) [23]. The results show that 
the WF of single photon-added squeezed thermal state 
(PASTS) always has negative values at the phase space 
center. The decoherence effect on the PASTS is exam- 
ined by the analytical expression of WF. It is found that 
a longer threshold value of decay time is included in sin- 
gle PASTS than in single-photon subtraction squeezed 
thermal state (STS). In this paper, as a natural exten- 
sion, we shall introduce a kind of nonclassical state — 
photon-addition two-mode STS (PA-TMSTS), generated 
by adding photons to each mode of two-mode STS (TM- 
STS) which can be considered as a generalized bipartite 
Gaussian state. Then we shall investigate the entangle- 
ment and nonclassical properties. 



This paper is organized as follows. In Sec. II we in- 
troduce the PA-TMSTS. By using the P-representation 
of density operator of thermal state, we derive the nor- 
mal ordering and anti- normal form of the TMSTS, which 
is convenient to obtain distribution function, such as Q- 
function and WF. Then a compact expression for the 
normalization factor of the PA-TMSTS, which is a Jacobi 
polynomial of squeezing parameter r and mean number 
n of thermal state. In Sec III, we present the nonclassical 
properties of the PA-TMSTS in terms of cross-correlation 
function, distribution of photon number, antibunching ef- 
fect and the negativity of its WF. It is shown that the WF 
lost its Gaussian property in phase space due to the pres- 
ence of two- variable Hermite polynomials and the WF of 
single PA-TMSTS always has its negative region at the 
center of phase space. Then, in Sees. IV and V are de- 
voted to discussing the entanglement properties of the 
PA-TMSTS by Shchukin-Vogel criteria and the quantum 
teleportation. The conclusions are involved in Sec. VI. 



II. PHOTON-ADDITION TWO-MODE 
SQUEEZED THERMAL STATE (PA-TMSTS) 

As Agarwal et al [2j . introduced the excitations on a 
coherent state by repeated apphcation of the photon cre- 
ation operator on the coherent state, we introduce theo- 
reticahy the photon-addition two-mode squeezed thermal 
state (PA-TMSTS). 

For two-mode case, the photon-added scheme can be 
presented by the mapping p -^ a^'^b^'^pa'^b^'^. Here we 
introduce the PA-TMSTS, which can be generated by 
repeatedly operating the photon creation operator a^ and 
b^ on a two- mode squeezed thermal state (TMSTS), so 
its density operator is 

/-^ = N-^ua^'^b^"S (r) pthiPth2S^ (r) a'^b", (1) 

where tti, n are the added photon number to each mode 
(non- negative integers), and Nm,n is the normalization 
of the PA-TMSTS to be determined by trp^"^ = 1, and 
S{r) = exp[r{a^b^ — ab)] is the two- mode squeezing oper- 
ator with squeezing parameter r. Here pthi,2 is a density 
operator of single-mode thermal state. 



n=0 V* "•" '■) 



(2) 



where n is the average photon number of thermal state 
Pthj {j = 1,2). For simplicity, we assume the average 
photon number of pthj {j = 1,2) to be identical. In 
addition, the P-representation of density operator pt^j 
can be expanded as [25| 



pthj 



n J TT 



(3) 



which is useful for later calculation and here \a) is the 
coherent state. 



A. Normal ordering and anti-normal form of the 
TMSTS 



In order to simplify our calculation, here we shall 
derive the normally ordering form of the TMSTS. For 
this purpose, we examine the two- mode squeezed co- 
herent states S \a^ P) (|q^,/3) = \a) |/3)). Note that 
\a) = exp[— I \a\ + aa^ |0) and the following transfor- 
mation relations Ea, S3 ' 



S{r)a^ S^{r) — a^ coshr 
S(T)b^S\T) = 6^ coshr - 



- b sinh r, 
a sinh r. 



(4) 



we see 



5'|a,/3) = sechrexp 



-i(|ar 



X exp \oL (a^ cosh r — b sinh r) ] 
X exp[/3 {\r cosh r — a sinh r) ] 



X 



exp [a'^^Uanhr] |00) , 



(5) 



where S{r) |00) =sechrexp [a'''6''' tanhr] |00) is used. 

Further noting e^^a^e~^^ = a^ ^ r^ and for operators 
A,B satisfying the conditions [A, [A, 5]] = \B, [A, 5]] = 
0, we have e^+^ = e^e^e-t^'^]/^ = e^e^et^'^]/^, thus 
Eq.© can be put into the following form 



5'|a,/3) = sechrexp 



-^(i«r 



^) — a/3tanhr 



X exp [(a'f'a + b^ P) sechr + a^b^ tanhr] |0((i^) 

Thus inserting Eq.(j6]) into Eq.(|3]) and using the vacuum 
projector |00) (00| =: exp[— a'^'a — 6'''6] : (where : : de- 
notes the normally ordering) as well as the IWOP tech- 
nique [28|, [29| , we can obtain 

P^ = SpthlPth2S^ 

= 4/"^^e-*(H^+l''l=)5|a,/?)(a,/?|5t 
n^ J TT^ 

= Al : exp [A2 {a'^b'^ + ab) - A3 {a'^a -|- b'^b)] : (7) 
where we have set 



Al 

A2 



sech^r 



(n + 1) — n^ tanh^ r 
(2n + 1) sinh r cosh r 

(2n + 1) cosh^ r -\-n'^^ 



cosh 



n cosh 2r 



(2n + 1) cosh^ r ^n^ 
and used the integration formula [30| 



(8) 



,C\z\'+e,z+r^z* 



-e-T',ReC<0. (9) 

Eq.© is just the normally ordering form of TMSTS to 
be used to realize our calculations below. 

In addition, using Eqs.©, (J9]) and the formula con- 
verting any single-mode operator O into its anti-normal 
ordering form [31|], 



O 



^■n<- 



-z\d\z)e 



\z\^-\-z* a—za^ -\-a^ a'. 



(10) 



where \z) is the coherent state, and the symbol : : de- 
notes antinormal ordering, (note that the order of Bose 

operators a and a^ within : : can be permuted), one can 
obtain the anti- normal ordering form of the TMSTS, 



p = Al'. exp 
where we have set 

Al = ■ 

A2 
As 



A2 {a'^b'^ + ab) - A3 {a'^a + b'^b)] :, (11) 



(n + 1)^ - (2n + 1) cosh^ r 

(2n + 1) sinh r cosh r 
(n + 1)^ - (2n + l)coshV 

sinh r -\- n cosh 2r 
(n + 1)^ - (2n + l)coshV 



(12) 



Ea. (fTT1l implies that the P function P(a,/3) of the TM- 
STSis 

P{a, 13) = Ai exp [ia (a*/3* + a/3) - Asi\af + |/3|^) 

(13) 
which leads to the P representation of density operator 
SpthiPth2S^ i.e., 



/ 



f^P{a,l3)\a,p){a,p\. (14) 



In particular, for the case without squeezing, r = 0, 
then Eqs.([7j) and (pT]) just reduce to, respectively. 



p' (r = 0) 



1 


: exp 


a'^a^b'^b' 
n + 1 


(n + lf 


1 . 
^:exp 


-4(« 


ta + 6t6) 


•7 



(15) 



as expected [25]. It is interesting to notice that, for the 
case of n = 0, corresponding to the two- mode squeezed 
vaccum state (TMSVS), Eqs.Q and ([II]) become 

/ (n = 0) 
= sech^r: exp [(^a^b^ + ab) tanhr — (^a^a + b^b)] : 

= —eschar: exp [a^a + b^b — (a^6^ + ab) coth r] :,(16) 

which are just the normal ordering form and anti- normal 
ordering form of the TMSVS. The second equation in 
Eq. (fT6|) seems a new result. Here, we should mention that 
the normal (anti-)normal ordering forms of the TMSTS 
are useful to higher-order squeezing and photon statistics 
MM for the TMSTS. 



B. Normalization of the PA-TMSTS 

To fully describe a quantum state, its normalization is 
usually necessary. Using Eq.([7j), the PA-TMSTS reads 
as 



SA 



Ai 



Nrr 



' /7t^?.t^p^2(a^6^+a6)-A3(a^a+6^6) mjLU . 



(17) 
Thus using the completeness relation of coherent state 
/ (fadPp |a, P) (a, /3| /tt^ = 1 and Eq.®, the normaliza- 
tion factor Nm,n is given by (Appendix A) 



Nrr 



where 



^ "^ "^ BliTt + T't'^+B^JTT' + tt') 



dT'^dt'^dT"^dt" 



(18) 



Here we introduce a new expression of generating func- 
tion for Jacobi polynomials in form (Proof see Appendix 
B) 



a 



2m+2n 



j^n-m ^^2 _ ^2^ 



^A{T't'+Tt)+B{TT'+t't) 



m\n\ < 



^(n— m,0) 



fB^±A^\ 



m ^ n 



(20) 



thus the normalization factor Nm,n can be put into (with- 
out loss of generality assuming m ^ n) 

Nra,n = mlnlB'^-'^LO^P^'"-^^ (^) , (21) 

where we have used the property of the Jacobi polyno- 
mials Pt'^\-x) = (-l)^Pi^'"^(x),and 

uj = fi^ -\- (2n + 1) cosh r, 

V = n (n + 1) cosh 4r + (n + cosh r)cosh2r. (22) 



Eg. ([21]) indicates that the normalization factor is related 
to the Jacobi polynomials, which is important for further 
studying analytically the statistical properties of the PA- 
TMSTS. Note Eq.([2I]) exhibits the exchanging symme- 
try. 

It is clear that, when tti = n = 0, Eq. ([2T]) just reduces 
to the TMSTS due to P^^'^^ {x) = I; while for n 7^ 
and m = 0, noticing Pq '^^ (x) = 1, Eq.([2T]) becomes 
^o,n = nlBi. For the case m = n, Nm,m is related 
to Legendre polynomial of the parameter - , because of 

Pn {x) = Pn{x)^ Pq{x) = 1. In addition, when n = 
leading to cj = Pi = cosh r and - = cosh2r, then 
Eq.([2l]) reads 



Nm^n (n = 0) = m\n\ cosh^^ rP^^'^-^^ (cosh 2r) , (23) 

which is just the normalization of two- mode photon- 
added squeezed vacuum state [3^]. 



III. NONCLASSICAL PROPERTIES OF THE 
PA-TMSTS 

In this section, we shall discuss the nonclassical proper- 
ties of the PA-TMSTS in terms of cross-correlation func- 
tion, photon statistics, anti-bunching effect and the neg- 
ativity of its WF. 



A. Cross-correlation function of the PA-TMSTS 



Pi = cosh r + n cosh 2r, 
P2 = (2n + 1) sinhr coshr. 



The cross-correlation between the two modes reflects 
(19) correlation between photons in two different modes. 
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FIG. 1: (Color online) Cross-correlation function between the 
two modes a and 6 as a function of r for different parameters 
(m,n) and n — 0.01. 



always larger than zero, thus there exist correlations be- 
tween the two modes. This implies that the nonclassical- 
ity is enhanced by adding photon to squeezed state. For 
given (tti, n) and n values, gm,n increases as r increasing; 
while gm,n decreases as n decreasing for a given {m^n) 
value. It is interesting to notice that for single-photon- 
addition TMSTS, the gm,n parameter presents its max- 
imum value, which implies that single-photon-addition 
TMSTS may possess a stronger nonclassicality than the 
other TMSTSs. To compare the further nonclassicality 
of quantum states for a different number added case, the 
measurments based on the volume of the negative part of 
the Wigner function [36], on the nonclassical depth [33|, 
and on the entanglement potential [38|, Vogel's noncalssi- 
cality criterion [39] and the Klyshko criterion [40] may be 
other alternative methods. 



which plays a key role in rendering many two-mode ra- 
diations nonclassically. From Eqs. (pT]) and (|2T]) we can 
easily calculate the average photon number in the PA- 
TMSTS, 



(^a^a) 



N. 



m-\-l,n 



N„ 



i,m 



N, 



m,n-\-l 



Nrr 



and 



{ah'^ab) 



-^^ m,.ii, 



(24) 



1. (25) 



Thus the cross-correlation function g^ n can be obtained 
by [25.] 



9mAr) 



1 



{a^bUb) 
(ata) (5t6) 



N, 



N„ 



m,n+l-^^m+l,n 



(iv, 



m,n+l 



iVm.n) (N, 



m+l,n 



■^^m^n) 



.(26) 



The positivity of the cross-correlation function Qm^n 
refers to correlations between the two modes. In par- 
ticular, when m = n = {) corresponding to the 
TMSTS, noticing 7Vo,o = l,A^o,i = ^"1,0 = ^i, 
and Ni^i = v^ then Eq.([26|) reduces to ^0,0 (^) = 
(2n + 1) sinh rcosh r/{Bi — l) , which implies that 
the parameter go^o{r) is always positive for any n and 
non-zero squeezing {Bi ^ 1). Further, for the case of 
^ = O5 ^0,0 (^) = coth^ r, which is just the correlation 
function of the TMSVS; while for r = 0, i.e., the TM- 
STS, ^0,0(0) = O7 so there is no correlation between two 
thermal states, as expected. On the other hand, when 
m = 0,n = 1, noticing Ni^2 = ^1 (3^; — a;) , 7Vo,2 = 
2Bl and P^^'^^ (x) = {3x - l)/2, then Eq.® be- 
comes go,i{r) = (v-u) Bi/[{2Bi - 1) (^; - ^i)]. Notic- 
ing that V — Bi > and {2Bi — 1) > 0, and v — uj = 
^ (2n + 1) sinh^ 2r ^ 0, so ^0,1 (^) is always positive. 

In order to see clearly the variation of ^^^^-parameter, 
we plot the graph of gm,n as the function of r for some 
different (m, n) and n values. It is shown that gm,n are 



Distribution of photon number of the 
PA-TMSTS 



In order to obtain the photon number distribution 
(PND) of the PA-TMSTS, we begin with evaluating the 
PND of TMSTS. For two- mode case described by den- 
sity operator p*^, the PND is defined by V{ma^ni) = 
{rria^nijl p^ \ma^rH)) . Employing the non- normalized co- 
herent state I a) = exp[aa'''] |0) leading to |n) = 

~T^\ ~h^ l^)| -0 {^1^ \^) ~ ^^^ )' ^^ ^^^^ ^^ ^^^ normal 
ordering form of p^ in Eq. ([7j) , the probability of finding 
(rria^nh) photons in the two-mode field is given by 



V{ma,nb) 



Ai 



^2ma+2nh 



mjnbl da'^-da^-dp^^dp*^^ 



X e 



(l-A3)(a*a+/3*/3)+A2(a/3+a*^*) 



= Ai[n(n + 1)]' 



p" 






where in the last step, we have used the new formula in 
Eq.(|20l), and 



u = (2n+l)cosh r-\-n'^^ 

p = (2n+l)cosh r — (n + 1) , 

((2n + 1) sinh 2r) V 4n2 (n + 1)^ 
^ ~ ((2n + l)sinh2r)^ -4n2(n + l)^ 



(28) 



Thus the PND of TMSTS is also related to Jacobi polyno- 
mials of the parameter x- In particular, when n ^ lead- 
ing to X — > 1, corresponding to the two- mode squeezed 



vacuum, Eq.(j27j) reduces to 

Vn^o{ma,nb) 



lim 



rrialnbl 



^^ [n2 + (2n + l)coshV] 



2 -irria+nb + l 



Ynin[ma,na 

E 

1=0 



(2n + 1)"' [n (n + i)]"»«+«^-2' ginh^' 2r 



22' (H)" (nfc - 0! (m„ - 0! 



tanh 



cosh r 



(29) 



which is just the PND of two- mode squeezed vacuum 
state [24]. On the other hand, when r ^ correspond- 
ing to the case of two-mode thermal state, leading to 
X -^ -1, i^ ^ (n + 1)^ , /i ^ -n^, Ai -^ 1/ (n + 1)^ and 



noting P^ j 



(nb-ma,0) 



(-1) = (-1)"^% thus Eq.(j27j) becomes 



Vr^oima^rib) 



(n + ir^^(n + l) 



rua + l ' 



(30) 



which is just the product of PNDs of two thermal fields, 
as expected. 

Using the result (|27l) and noticing a'^b'^ \ma,nb) = 
^^ma\nb\/{ma - m)\{nb - n)\ \ma - m.rib - n), we 

can directly obtain the PND V^^{ma^ni)) = 
{ma,nt\p^^\ma,nb) of the PA-TMSTS as 

N-l^mglribl 
{rria — m)\{ni) — n)\ 
X {rria — m^ni) — n\ p^ \ma — m^nij — n) 



{rua -m)\(nh - n)\ 



V{ma -m.nij -n). (31) 



Eq. (|3T]) is a Jacobi polynomial with a condition rria ^ 
m and niy ^ n which shows that the photon-number 
{rria^ni)) involved in PA-TMSTS are always no-less than 
the photon-number {m^n) operated on the TMSTS, and 
there is no photon distribution when rria < ^ and n^ < 
n. Here we should point out that this result (|27|) can 
be applied directly to calculate the PND of some other 
non-Gaussian states generated by subtracting photons 
from (or adding photons to) two-mode squeezed thermal 
states, such as a'^b'^p^a^'^b^'^, and a'^b^'^p^b'^a^'^. 



C. Antibunching effect of the PA-TMSTS 

Next we will discuss the antibunching for the PA- 
TMSTS. The criterion for the existence of antibunching 
in two-mode radiation is given by [41j 

n.^'''rtr''^-^<o- (32) 



In a similar way to Eq. (|26|) we have 

-1-^171,71 -'-^171,71 

(5t2;,2^ = i^,n+2_4i^,n+1^2. (33) 

Thus, for the state p*^^, substituting Eqs.(j33]), (J25j) and 
(f2T]) into Eq.(j32j), yields 



A/'^+2,n + A/'^^n+2 + 2 (r^ — 7Vm+l,n+l) /o/,\ 
^ab = 57TT — TTT d^V 

(^ = Nm,n — Nm+l,n " Nm,n+l) • 

In particular, when m = n = (corresponding to the 
TMSTS) leading to Nq^q = l,A^o,i = A^i,o = ^i, and 
Ni^i = V, 7Vo,2 = A^2,o = Bf, thus Eq.(j34]) becomes 



R 



'ab,m=n- 



(2n + 1) (4 cosh 2r + (2n + 1) sinh^ 2r) 

~ (2n + 1) [(2n + 1) cosh 4r - 2 cosh 2r] + 1 * 

(35) 

From Eq.(j35j), it is easily seen that Rab^m=n=o < for 
any n and non-zero r values. In addition, when m = n, 
the PA-TMSTS can always be antibunching for a small 
value n (see Fig. 2 (a)). However, for any parameter values 
m^n{m 7^ n), the case is not true. The Rab parameter 
as a function of r and m^n is plotted in Fig. 2. It is 
easy to see that, for a given m the PA-TMSTS presents 
the antibunching effect when the squeezing parameter r 
exceeds to a certain threshold value. For instance, when 
m = and n = 1 then Rab may be less than zero with 
r > 0.1 thereabout (n = 0.01). The value Rab parameter 
increases with n increasing. 



D. Wigner function of PA-TMSTS 

For further discussing the nonclassicality of PA- 
TMSTS, we examine its Wigner function (WF) whose 
partial negativity implies the highly nonclassical prop- 
erties of quantum states. In this section, we derive the 
analytical expression of WF for the PA-TMSTS. The nor- 
mally ordering form of the PA-TMSTS shall be used to 
realize our purpose. 

For the two-mode case, the WF W (a, /3) associated 
with a quantum state p can be derived as follows [42|: 



W 



X exp [2 {azl - a*zi) + 2 {f3z^ - /3*Z2)] (36) 



where \zi^Z2) = \zi) \z2) is the two-mode coherent state. 
Substituting Eq. (fT7j) into Eq. (|36|) . we can finally obtain 
the WF of the PA-TMSTS (see Appendix C), 



2(ata6t6) 



Wm,n {a, P) = Wo (a, P) Fm,n {a, P) , 



(37) 




FIG. 2: (Color online) Rat as a function of r for different 
parameters (m,n) and n = 0.01. 



where V^o (o^, P) is the WF of TMSTS, 

"^ ^ _cosh2r 



Wo {a, 13) = 
and 

Fm,n {a, P) 



(2n + 1) 
„ sinh 2r 



exp 



2n+l 



2n+l 

{Pa + a* 13*) 



(38) 






H, 



m—l,n—j 



Ri Rs 



i\fK^ iJKk 



(39) 



where we have set 
Ri 



n + cosh r sinh r cosh r 

^1 = — ^_ , , — ,^3 = 



2n + l 



2n + l 



(40) 



Equation (|37|) is just the analytical expression of the WF 
for the PA-TMSTS, a real function as expected. It is 
obvious that the WF lost its Gaussian property in phase 
space due to the presence of two- variable Hermite poly- 
nomials Hm-l,n-j {X,y). 

From Eq. (|39|) , we see that when m = n = corre- 
sponding to the TMSTS, Fo,o = 1, and Wm,n{(^,P) = 
Wo (a,/3); whereas for the case of m = and n ^ 0, 




FIG. 3: (Color online) The Wigner function W(a, /3) in phase 
space ((3+,P+) for several different parameter values (?7i,n) 
and n with r = 0.3. (a) m=0,n=l,n = 0.2; (b) m=0,n=l,n = 
1; (c) m=n=l,n = 0.2 and (d) m=n=l,n = 1. 



noticing Ho,n (^, y) 
duces to 



y"" and A^o,n = n\B]^, Eq.(|37D re- 



Wo,n {a, p) = Wo (a, f3) {-K^/B^f i^ndi^sl' /i^i), 

(41) 
where Ln is the n-order Laguerre polynomial. Eq. (|^T]) is 
just the WF of the PA-TMSTS generated by single-mode 
photon addition, which becomes the WF of the negative 
binomial state S{r) |0,n) with n = [JOSAB, ??]. In 
particular, for the case of single photon-addition, n = 1, 
it is found that 



^0,1 (0,0) 



(n + cosh^ r) / (2n + 1)' 
(cosh 



ncosh2r) tt^ 



(42) 



which implies that the WF of single PA-TMSTS always 
has its negative region at the phase space center a = 
P = 0. The maximum value of |I^o,i| decreases with 
the increasement of n and r but not disappears, which 
can be seen clearly from Fig3,4. (a) and (b). Further, 
there are more visible negative region than the WF for 
the case of tti = n = 1. And the negative region will 
be absence for the latter with the increasing n value (see 
Fig3,4. (c) and (d)). In addition, from Figs 3 and 4, the 
squeezing in one of quadratures is clear, which can be 
seen as an evidence of nonclassicality of the state. For a 
given value m and several different values n (7^ m), the 
WF distributions are presented in Fig. 5, from which it is 
interesting to notice that there are around \m — n\ wave 
valleys and \m — n\ -\- 1 wave peaks. 



IV. ENTANGLEMENT PROPERTIES OF THE 
PA-TMSTS 

It is well known that photon subtraction/addition can 
be applied to improve entanglement between Gaussian 
states 115|,|43|, loophole- free tests of Bell's inequality [4^], 




FIG. 4: (Color online) The Wigner function W(a, /3) in phase 
space (Q-,P-) for several different parameter values (m,n) 
and n with r = 0.3. (a) m=0,n=l,n = 0.2; (b) m=0,n=l,n = 
1; (c) m=n=l,n = 0.2 and (d) m=n=l,n = 1. 




FIG. 5: (Color online) The Wigner function W(a, /3) in phase 
space ((5+,P+) for several different parameter values (m,n) 
with n = 0.2 and r = 0.3. (a) m=0,n=2; (b) m=l,n=2; (c) 
m=l,n=3, and (d) m=2,n=3. 



and quantum computing [19]. In this section, we exam- 
ine the entanglement properties of PA-TMSTS only with 
single and two photon-addition. Here, for a bipartite con- 
tinuous variable state, we shall take the Shchukin-Vogel 
(SV) [39| criteria to describe the inseparability of PA- 
TMSTS. 

According to the SV criteria, the sufficient condition 
of inseparability is 



SVm,n = (a^a-^)(b^b-^ 



(a^6^)(a6) < 0. (43) 



where we have set 



t,T,t'T'=0 






5! (g — r + s)\ (r — 5)! (/ + r 
Thus SV is given by 



sy. 



N 2 \ N 



-'-^m,n 



-1-^171,71 



(46) 



Next, we examine two special cases. For the case of ?7i = 
0,n = 1, using Eqs. (|2T]) and (|^5]) , as well as noticing 
A^o,i,i,2 = ^^1,0,2,1 = 251^2, Eq.(j46l) becomes 



^^0,1 



V 

B^ 



2Bi 



ABl 



(47) 



While for the case of tti = n = 1, it is shown that 
(iVi,2,i,2 - iV2,i,2,i = 2 {2Bl + Bfi B2) 



.^,.(s,(3-^)-|) = -4 



(2i?? + B2)2^2 



■ (48) 



In particular, when n = 0, i.e., the single PA-TMSVS, 
Eq. (|T7|) is always negative for any r > 0, as expected 
(also see Fig. 6 (a)). In general, it is difficult to obtain 
the explicit expressions of the sufficient condition of in- 
separability for the above cases. Here, we appeal to the 
number calculation shown in Fig. 6. It is shown that for 
single PA-TMSTS with a smaller average photon num- 
ber n, the condition SVo^i < can always be satisfied 
only if r > 0; while for a larger n then the condition 
SVo^i < is satisfied only when the squeezing parame- 
ter r exceeds a certain threshold value r^. However, it is 
very interesting to notice that for the photon-subtraction 
TMSTS, there is a threshold value Tc for any n, i.e., 
r > re = ^ ln(2n + 1) [45], which is different from the 
case of single PA-TMSTS. For instance, for n = 1, the 
two threshold values are r^ ~ 0.31 and Tc ~ 0.55. This 
comparision may imply that the photon-addition to the 
TMSTS can be more effective for the entanglement en- 
hancement than the photon-subtraction from the TM- 
STS. On the other hand, for the case of the PA-TMSTS 
with m = n = 1 (see Fig. 6 (b)), it is found that a 
certain threshold is needed for satisfying this condition 
SVi^i < 0, which is also smaller than that of the photon- 
subtraction TMSTS. 



In a similar way to derive the normalization factor 
Eq.dni, using Eqs.([I71) and ([IH]), we have 



QUANTUM TELEPORTATION WITH 
PA-TMSTS 



(at^t) 



N, 



m,m+l,n,n+l 

~N ' 

-'■^ m,.n 



(ab) 



N^+ 



l,m,n+l,n 



A^^ 



, (44) 



As mentioned above, photon-subtraction from or 
photon-addition to bipartite Gaussian states can be used 




FIG. 6: (Color online) The sufficient condition of insepara- 
bility as the function of r and n, for (a) m=0 and n=l; (b) 
ni=n=l. 



to improve the entanglement. In this section, we investi- 
gate the quantum teleportation with PA-TMSTS, espe- 
cially for the cases m = 0,n = 1 and m = n = 1. The 
role of teleportation in the CV quantum information is 
analyzed in the review Ref.[46]. 

Here, we consider the QT by using PA-TMSTS as 
entangled resource. Using the normal ordering form 
Eq.([T7j) and noticing the displacement operator Da (a) = 
g|a| /2^-a*a^aa ^ ^^^^1 the characteristic function (CF) of 
PA-TMSTS is given by (see Appendix D) 



XE{a,f3) 



-(Bi-i)(|a|2 + |/3|2)+B2(a/3+a*/3*) 



Nrr 



X- 



d' 



i2m+2n 



^gt(aBi-B2/3*)+r(/3B2-Bia*) 
^ ^T'{aB2-Bi(3*)^t'{(3Bi-B2a*) 



^Bx{T't'+tT)+B2{t'tWr) 



t,T,t'T' = 



m 



where for further calculation the differential form of xe 
is kept. 

To quantify the performance of a QT protocol, 
the fidelity of QT is commonly used as a measure, 
T =tT {pin pout) 1 a overlap between a pure input state 
pin and the output (teleported, mixed) state pout- For 
a CV system, a teleportation protocol has been given in 
terms of the CFs of the quantum states involved (input, 
source and teleported (output) states) [47]. It is shown 
that the CF Xout {v) of ^he output state has a remarkably 
factorized form 



Xoutiv) =Xin{v)XE{v'',V)' 



(50) 



where Xin{v) ^^^ XEiv'^^v) ^^^ the CFs of the input 
state and the entangled source, respectively. Then the 
fidelity of QT of CVs reads [47|] 



T 



(Pt] 



Xin {r]) Xout {-r]) • 



(51) 



Here, we consider the Braunstein and Kimble proto- 
col ^] of QT for single- mode coherent-input states I7). 
Note that the fidelity is independent of amplitude of the 
coherent state, thus for simplicity we take 7 = 0, then 
we have only to calculate the fidelity of the vacuum input 
state with the CF Xin iv) = exp[— I77I /2]. On substitut- 
ing these CFs into Eq. ([5T]) , we worked out the fidelity for 



teleporting a coherent state by using the PA-TMSTS as 
an entangled resource. 



*/ rr 



(m + n)! (Bi + B^) 



'm-\-n 



[(2f 



Bo 2™+"+iA^„ 



l)e 



2r 



1] 



m^n 



(m + n)\ 



(2n + 1) e-2^ + 1 22^+2^7V^,n 



. (52) 



It can be seen that the fidelity is not only dependent on 
the parameter r, the average photon-number n, but also 
on the photon number (772, n) added to each mode of the 
TMSTS. In particular, when m = n = 0, Eq.(j52j) just 
reduces to 



1 



^'° (2n + l)e-2-. 



(53) 



which leads to the condition r > ^ In (2n + 1) for satis- 
fying the effective QT with T > \ which is the classical 



limit. In addition, for the case of n 
added TMSVS, Eq.(j52]) becomes 



0, i.e., the photon- 



x-o _ 



(e^-- + 1) 



n-\-m 



{m -\- n)\ 



-2r ^ I 22^+2^ A^^ 



(54) 



Further, when {m,n) = (0,0), (1,1) and (0,1), Eq.(j53l) 
just reduce, respectively, to 



-^0,0 

-^1,1 

-^0,1 



(l + tanhr)/2, 

(1 + tanhr)^ 

4(1 + tanh^ r) ' 

1 + tanh r 



4(1 - tanhr) 



sech^r. 



(55) 



The two expressions Tq q and Ti ^ are agreement with 
Eqs.(15) and (17) in Ref. [49]. 

In Fig. 7, for some given {m^n) values, the fidelity of 
teleporting the coherent state is shown as a function of r 
by using the PA-TMSTS as the entangled resource. It is 
shown that the fidelity with this resource is smaller than 
that with TMSTS, although the PA-TMSTS posesses 
larger entanglement [49[. In addition, for the symmet- 
rical case m = n^ when the squeezing parameter r ex- 
ceeds a certain threshold value, the fidelity increases with 
a increasing m (see Fig. 7(c)); while for non-symmetric 
case m ^ n, the fidelity decreases with increasing n (see 
Fig. 7(b)). For the former, the threshold value r decreases 
with increasing m{= n); the case is not true for the lat- 
ter. This indicate that the symmetrical PA-TMSTS may 
be more effective for QT than the non-symmetric case. 



VI. CONCLUSIONS 

In this paper, we introduce the PA-TMSTS and inves- 
tigate its entanglement and nonclassicality. By using the 
coherent state representation of thermal state, the nor- 
mally and antinormally ordering forms of the TMSTS 





FIG. 7: (Color online) The fidelity as the function of r for 
several different (rQ,n) values and n= 0.01. 



are obtained. Based on this, the normalization factor of 
the PA-TMSTS is derived, which is related to the Jacobi 
polynomials of the squeezing parameter r and average 
photon number n of the thermal state. Then we dis- 
cuss the nonclassical properties by using cross-correlation 
function, distribution of photon number, antibunching ef- 
fect and the negativity of its WF. It is found that the WF 
lost its Gaussian property in phase space due to the pres- 
ence of two- variable Hermite polynomials and the WF of 
single PA-TMSTS always has its negative region at the 
center of phase space. Further, there are more visible 
negative region than the WF for the case of m = n = 1. 
And the negative region will be absence for the latter with 
the increasing n value. The entanglement properties of 
the PA-TMSTS by Shchukin-Vogel criteria and the quan- 
tum teleportation. It is shown that the photon-addtion 
to the TMSTS can be more effective for the entangle- 
ment enhancement than the photon-subtraction from the 
TMSTS; And using the PA-TMSTS as an entangled re- 
source, the fidelity for teleporting a coherent state is not 
only dependent on the parameter r, the average photon- 
number n, but also on the photon number (tti, n) added 
to each mode of the TMSTS. From this point, the sym- 
metrical PA-TMSTS may be more effective for quantum 
teleportation than the non-symmetric case. 
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Appendix A: Derivation of Eq. ([T8|) 

According to the normalization condition, trp*^^ = 1, 
we have 



j(fa(fp\a,p) (a,/3|/7r2 = 1 and Eq.®, Eq.(Al) 



N, 



-ad P |^|2m|^|2n^A2(a*^*+a/: 



dT'^dt'^dT"^dt" 



I 



e 

(Pa(Pf3 



Asi\a\^ + \i3\^) 



■ exp 



-As{\a\ 



') 



+ {A2I3 + T)a+ {A2I3* +t)a*+ t'i3 + t' (3* ,^ ^ ^, ^,^^ . 

'(A2) 



Using Eq.©, (A2) becomes 



N, 



Q^mQfmg^ingfn ^^ 



£l3 



■ exp 



13- 



'A2t 

= LHSofEq.(Il8l 



^ + Or + 



r Ai 


-^,«^ 


[ Az "^' 


Tt' 


t,T,t',T'^0 




(A3) 



where we have used Ai/{A\ — A2) = 1 and B2 = A2 



/{Al- 
AzIiAl 



Al) 



Al). 



(2n + 1) sinhr coshr as well as B\ 



Appendix B: New expression of generating 
function for Jacobi polynomials 

In this appendix, we shall prove Eq. (|2Q|) . Rewriting 



R 



^2m+2n 



^A(t' t' ^Tt)^B{TT' ^t' t) 



Expanding the exponential items, we see 

^^+J^s+fc g2m-\-2n 



(Bl) 






min[m,nj 

E 

s=0 



(m!n!)^A^+^ 



^k+jp+s^n+k^n+s 



t,T,t',T'^Q 



s\s\ {n — s)\ {m — s)\ 



A^ 



(B2) 



Comparing Eq.(B2) with the standard expression of Ja- 
cobi polynomials. 



Pi".«(x) 



1 



E 

/c=0 



m-\- a 
k 



1 



we can find that taking m 



H 



Nm,n = AitT 



imi-\nA2(^a'^b'^-\-ab^—A3(^a'^a-\-b'^b^min . 



A i^2m-\-2n 

dT'^dt'^dT''^dt' 



-tr 



E 

s=0 

m\n\ 
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mir[[rn, 

(Al) X Y. 



(m!n!)^A^+^ 
s\s\ {n — s)\ {m — s 

-J \ —m 

y \ Am-\-n 

i\n\ 



A^ 



X — 1 

(B3) 
n and y = 



{(^ 



Using the completeness relation of coherent state 



fe=0 

= m\n\A" 



k\k\ {n - ky. {m - ky. \y - 1 



(B4) 
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In a similar way, for n ^ tti, we also have 

H = mlnlA'^-'' {B^ - A^)"" P^^-^'^^ {y) . (B5) 

Thus we finish the proof of Eq. (|2Q|) . 

In addition, when m = n^ Eq. (|2Q|) becomes 



/_-| \m-\-n 



:^2m+2n 



K- 






1^3 



d^ 






g^mQfmg^fmgfm 



^A(T't'+Tt)+B(TT'+t't) 

(^At' -\-Bt)t' -\-(^AT-\-Bt'y 



f^^\l U'l'y ^Rit+R^t' -RIt-RIt' +K:i{TT' +tt') 
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(-1) 
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i^i 



i+j 



Q21+2J 



^™," i^Q ^-J- d{-Rl)^ dR\d{-Riy dRi 
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d 
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T,t'=0 



= (m!)^(B2-A2)"P, 



B^ + A" 



(B6) 



where Pm (x) is the Tnth Legendre polynomials. Eq.(B6) 
is just a new formula for the generating function of Leg- 
endre polynomials P^(x), which is different from the new 
form found in Ref. [50] . In fact, one can check Eq. (B6) by 
expanding directly the whole exponential items and com- 
paring with the standard expression of Legendre polyno- 
mials. 

Appendix C: Derivation of Eq. (|37p 
Substituting Eq.i^T^ into Eq.(|36l) and usiing Eq.®, 
we have 



t,t',r,r'=0 

(C5) 

Further using the generating function of two-variable 
Hermite polynomials, 

fim fin 

^ ^ ^-Atv+Bt+Cv 



Q^m Qyn 



Eq.(C5) can be put into the following form 
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where Wo {p^^P) is defined in Eq. (|38|) . and 
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thus we can obtain Eq. (|39l) . 
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Appendix D: Derivation of Eq. (|49p 

Using the displacement operator Da (a) = 

well as the normally ordering form of PA-TMSTS (fTT]) . 
the CF of PA-TMSTS is given by 

XE{a,/3)=tr[Da{a)D,{/3)p'^] 



and 



Nrr 



Sa^a/3^a(-a*)'^S(-/3*)' 



Rl = 2{Kia-K3f3^) 
Rs = 2{Kip-Ksa'') 



-^2 5 



X tr 
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as well as 
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sinh r cosh ? 
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(C4) 



(Dl) 

In a similar way to derive Eq. (p!8|) . using Eqs.(Al) and 
(p!8|) . one can directly obtain 



Xe{(^,P) = 



.(|/3|2 + |a|2)/2 



d' 



\2m-\-2n 



Nm,n da^dp^'d (-a*)^" d (-/3*)^' 



Expanding the partial exponential items in Eq.(C2), then 
Eq.(C2) becomes 



X e 
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(D2) 
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Taking the following transformations 

a ^ a -\- r^ —a* -^ t — a* ^ 

which leads to 

^Bl(a(-a*)+^(-^*))+B2(a^+(-a*)(-^*)) 



(D3) 



-^ exp 



-B^{\a\' + \p\') + B2{ap + a*n 



X exp [Bi {r't' + tr) + B2 {ft + t't)] 

X exp [t {aBi - Ba/?*) + t {13 B2 - Bia*)] 

X exp [t' {aB2 - Bil3*) + t' {/3Bi - Baa*)] , 

thus Eq.(D2) becomes Eq. (|39|> . 



(D4) 
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